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Abstract. Electron-impact excitation of Xe atoms in pure Xe gas and in a 
Xe(10 %)— Ar(90 %) mixture has led to the discovery of infrared (IR) luminescence 
of Xe2 excimers. The investigation of the emission spectrum at low gas density has 
allowed the identification of the molecular states involved in the transition. When 
the gas density is increased to values up to 40 times larger than the density of the 
ideal gas at standard temperature and pressure, the interaction of the excimer 
with the dense environment produces a strong red-shift of the spectrum that is 
interpreted in terms of many— body effects. 
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1. Introduction 

The luminescence of rare gas excimers is studied extensively for applications in the 
vacuum ultraviolet (VUV) including production of coherent and incoherent VUV light 
sources [l| and detection of ionizing radiation Q . Xe is of particular interest because 
of its efficiency to convert the excitation energy in intense VUV radiation. When 
the gas is excited by means of several techniques including electrical discharges [U, 
irradiation with ionizing particles 0; HI; @| , or resonance lines 0; 13 E3 11 | , three- 



body collisions of excited and ground state atoms lead to the formation of the excited 
dimer, i. e., excimer, Xe£. The decay of the lowest lying excited levels towards the 
dissociative ground state of the molecule gives origin to the emission of the 1st and 
2nd VUV continua at about 152 and 170 nm, respectively. 

Owing to their relevance for so many applications, these VUV excimer bands 
have been thouroughly studied theoretically [HI; E3 0; EH and experimentally 
HI; E3 E3; 0; El in order to ascertain the molecular structure and potential energy 



curves and the kinetics of the deexcitation processes leading to the population of the 6s 
excited atomic states that are eventually involved in the production of these molecular 
states 0; B E3 IH [H; S3 H S3 S3 123 H- 



At quite low gas pressure (P < 20 kPa), the VUV luminesce signal consists of 
the 1st continuum. On the contrary, at higher pressure (P > 50 kPa), only the 2nd 
continuum is present in the spectrum. This different behavior of the two continua as a 
function of pressure has been explained by time-resolved spectroscopic studies showing 
that the 1st continuum is due to radiative transitions from the vibrationally excited 
(0j)i/>o state correlated with the resonant 5p 5 6s( 3 Pi) atomic state towards the 
dissociative 0+ ground state, whereas the 2nd continuum consists of the overlapping 
bound-free emission from the lowest vibrationally relaxed (0~,l u ) states correlated 
with the metastable 5p 5 6s ( 3 i"2) atomic state. 

Nearly no attention at all has been paid to the possibility that molecular 
infrared (IR) transitions may occur in the cascade of processes leading to the 
formation of the excimer states radiating in the VUV. Only an IR spectrum centered 
about 800 nm has been observed and its appearance has been attributed to a 
0+ (6p[l/2]o) — > 0+ (6s[3/2]i) 1 ,/ > o transition [9(. One possible reason for this lack 
of interest on IR processes may be that the minimum of higher-lying excimer states 
occurs at an internuclear distance at which the weakly bound molecular ground state 
is strongly repulsive and is therefore not so easily accessible by multiphoton selective 
excitation. On the contrary, broad-band excitation techniques using high-energy 



charged particles [6j; l2l| may produce excited atoms of such high kinetic energy that 
may collide with ground state atoms at short distance yielding molecular states in 
higher energy levels although their parity cannot be controlled. 

We have recently discovered a broad band in the near IR region centered at a 
wave number v m ~ 7800 cm -1 , i. e., A ~ 1.3 /im. The emission is produced by exciting 
the gas at room temperature by means of short pulses of electrons of 70 keV energy 
[29I ; I30I ]. Technical details of the experiment can be found in literature. 

The IR band is quite broad. At P w 20 kPa, its FWHM is T « 900 cm" 1 . The 
relative width T/P m m 0.115 is comparable with that of the 2nd continuum [3 l3lj ]. 
This IR band occurs both in the pure gas and in an Ar (90 %)-Xe (10 %) mixture. 
As the features of the spectrum in both gases in the limit of low P are identical, we 
attribute the emission to a Xe2 molecular bound-free transition. 

Improvements in our experimental apparatus have allowed us to carry out more 
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accurate time-integrated measurements of the excimer IR emission at low pressure 
and also to extend the investigated density range up to density nearly 40 times larger 
than the density of the ideal gas at standard temperature and pressure (273 K and 0.1 
MPa). By so doing, we have now the opportunity to better assign the molecular states 
involved in the transition and to elucidate the interaction of the excited molecules with 
a dense gas environment. 

The electronic structure of homonuclear rare gas excimers can be described by 
an ionic molecular core surrounded by an electron in a diffuse Rydberg orbital much 
larger in diameter than the internuclear distance [13; [H; [33| . Such a state can still 
exist in a high-pressure environment provided that the optically active electron is only 
weakly scattered off the atoms of the host gas. This occurs if the electron mean free 
path is larger than the orbit of the Rydberg state as it happens to be at low density. 
It is therefore interesting to look at the evolution of a molecular Rydberg state when 
the environment becomes increasingly denser. In our previous experiments [29H30T| we 
have mainly detected a red-shift of the IR emission spectrum of the excimer. In this 
experiment we confirm the previous results in an extended density range. 

2. Low density results 

In Fig. [1] we show an experimental IR emission spectrum recorded in pure Xe gas 
at P w 0.1 MPa at room temperature in an extended wave number range with a 
resolution of 16 cm -1 . The broad line shape is the molecular excimer band centered 
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Figure 1. High-resolution excimer IR emission spectrum. The atomic line 
related to a 7s — 6p transition is shown. 



at v m « 7800 cm" 1 . 

Several atomic lines appear in the emission spectrum as to be expected because of 
the broad band excitation induced by electron impact. It is particularly interesting to 
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note that an atomic line produced in a 7s [1/2]° — 6p[l/2] atomic transition [34[ 
appears to be embedded in the excimer band. The closeness of such an atomic 
transition to the excimer spectrum suggests that the upper molecular state may be 
correlated with a 5p 5 6p atomic limit. 

Only recently theoretical potential e nergy curves for the higher lying molecular 
states of Xe2 have appeared in literature 1J; LL5J . It is therefore now possible to test 
the accuracy of some of the theoretical potentials by comparing this first experimental 
excimer IR emission spectrum with that simulated by using standard Franck-Condon 
calculations [351: 

We thus assume that a bound-free transition is responsible of the observed 
emission and that the upper bound state is correlated with the 6p [l/2]o atomic state. 
The molecular potential correlated with this state shows an avoided crossing with 
the g molecular state related to the 5d[l/2]i state leading to predissociation [iij ]. 
The predissociation lifetime is estimated to be t v w 10~ 10 s, whereas the mean time 
between collisions at the density of the present experiment, N e w 2.4 • 10 25 m -3 , is 
i~c < 10~ n s. Thus, the large collision rate electronically stabilizes excimers before 
they predissociate. 

At present, there are no estimates of the radiative lifetime and of the decay rate 
for vibrational relaxation for this bound state. We thus assume that they are not too 
different from those of the states responsible for the emission of the VUV continua. 
The radiative lifetimes of the highly excited vibrational states of the 0+ and (1,,,, 0~) 
molecular states correlated with the 6s atomic limit are estimated to be t\ « 5 ns [201 ] 
and T2 ~ 40 ns [ld| . respectively. Similarly, the decay rate for vibrational relaxation 
k 3 is the same for all those state [§], yielding a decay time r 3 = (k^N^ 1 w 0.65 ns. 

This hierarchy of characteristic times r c < r p < T3 < t\ < T2 supports 
the assumption that collisions stabilize excimers electronically and quickly establish 
thermal equilibrium among their internal rovibrational degrees of freedom. 

Franck-Condon calculations require the knowledge of the potential energy curves 
whose choice among those reported in literature [151 ] has to comply with several 
requirements: 1) the upper bound state must belong to the 6p atomic manifold; 2) 
the parity of the initial and final state must obey the selection rules u <-> g and + ^ + ; 
3) the energy difference between the two potential energy curves at the coordinate of 
the minimum of the bound state potential must approximately be equal to the central 
wave number v m of the observed spectrum. 

According to these criteria, the upper bound state has been assigned as the (3)0+ 
ungerade state correlated with the 6p( 1 I?2) limit (with some contributions due to the 
mixing with a bd atomic state) The lower dissociative state has been assigned 

as the (1)0+ gerade state correlated with the 6s ( 3 Pi) resonant atomic state [15j . For 
both of these molecular states the projection Q, of the electronic angular momentum 
along the molecular axis is zero. Thus, as in this transition O = 0^O = 0, the 
additional selection rule for the total molecular angular momentum A J = ±1 must 
be also taken into account. 

The intensity of the line shape, at infinite resolution, is computed within the 
centroid approximation because the electronic transition moments are not known as 
a function of the internuclear coordinates for the states at hand thus yielding 



1(9) oc ^e-^^'i> 4 j(J' + l)|<e", J' + l\v' , J'} 

v'J' ^ 



+ 
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j'-l\v',J') (1) 



+ J' 

in which the selection rule J' — J" = ±1 is enforced. As usual, a single prime refers 
to the initial state (the bound one, in this case) and a double prime refers to the final 
state [33| . E v iji is the energy of a rovibrational level of the bound state and is given 

by 

E v ,j, = D' e [1 + € V ,} + B' e J' (J' + 1) (2) 

where D' e is the dissociation energy of the upper state and e v > are the scaled vibrational 
eigenvalues measured from the dissociation limit. B' e = h 2 /2m r R e > is the rotational 
constant, in which m r is the Xe2 reduced mass and R e i is the equilibrium internuclear 
distance of the bound state. The energy of the emitted photon is 

v = [T' e - T 1 : D'i + E v ,j, e"] (3) 

in which T' e and T" are the values of the minimum of the potential energy curves for 
the upper and lower states and D' e and -D" are their well depth. It has to be recalled 
that the 0+ state, though mainly dissociative, has a small attractive well at quite large 
internuclear distance, with the minimum located in R e n . 

\v' , J'} is the eigenfunction of a rovibrational state of the bound potential. |e", J") 
is a scattering state of kinetic energy e" and total angular momentum J" in the 
vibrational continuum of the dissociative potential. 

The Boltzmann factor in Eq. [T] accounts for the equilibrium thermal distribution 
of the rovibrational degrees of freedom with ~ 208.5 cm -1 at room temperature. 

The vibrational eigenvalues e v > and eigenfunctions \v', J') are found by integrating 
numerically the Schrodinger equation for the rotationless potential of the bound 
state by adopting the Numerov-Cooley finite difference scheme [37J and replacing 
the centrifugal potential with the constant value B' e J'(J' + 1) [38|. For numerical 
purposes, the literature 0+ potential has been accurately fitted to a Morse potential 

V +(R) =T' e + D' e {l - e [-^'( fl "^')l} 2 (4) 

with T' e = 13860 cm -1 , D' e = 1717 cm ~\ R e , = 3.23 A, and f3 e ,R e , = 6.734. The 
bound state accomodates approximately 34 vibrational states though only the first 
few contribute significantly to the spectrum because of the Boltzmann factor. 

The numerical value of the rotational constant is quite small B' e « 2.47 • 
10 -2 cm -1 , corresponding to a rotational temperature 9 r « 3.5 • 10~ 2 K. Thus, 
the thermal population of rotational states is non negligible up to J' « 250. 

The scattering eigenfunctions of the dissociative potential \e",J") are found by 
numerically integrating the Schrodinger equation with the effective potential 

Vf'VoiW + ^J'^ + l) (5) 

in which V + (R) is the rotationless potential of the dissociative state. V + (R) 

is characterized by a very shallow well of depth D" m 218 cm _1 occurring at 
R e n « 4.92 A and by a value T" w 4779 cm 1 15]. For numerical purposes, the 
potential is cast in the form V +(R) = T'J + D"f(x), in which f{x) is an analytical 
HFD-B potential [3^] and x — R/R e ii. A Runge-Kutta 4-th order scheme with 
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adaptive stepsize control has been adopted to integrate the Schrodinger equation [40j • 
The scattering wave functions are normalized to unitary incoming flux 41] 

^ = (R\e",J") [^j'\iu{kR + r 1 ) (6) 

in which k = [(2m r /S 2 )e"]^ 2 is the wave vector and r\ is the appropriate phaseshift. 

The overlap integrals in Eq. [1] are evaluated by spline interpolation and 
quadrature [13] . Moreover, the calculated spectrum has been convoluted with the 
instrumental function to account for the finite resolution of 16 cm -1 . In Fig. [2 we 
compare the experimental and calculated spectra. 
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Figure 2. Comparison of the experimental and calculated spectra. Dashed 
curve: spectrum computed by using literature potentials. Thick solid line: 
spectrum computed by using a i?~ 12 repulsive potential with adjustable 
parameters. 



The dashed line is the spectrum calculated by using the theoretical potential 
energy curves found in literature. The agreement is only qualitative. The synthetic 
spectrum is centered at much smaller wavenumbers than the experimental one and 
also its width is much larger than observed. These features are the consequence of 
the fact that the theoretical potential energy curves are too close to each other and 
that the dissociative potential is probably too steep in the region of the minimum of 
the bound state. The only correct feature shown by the computed spectrum is its 
asymmetric shape with a longer tail towards higher energies. This feature is due to 
the non negligible contribution of vibrational states with v' > 0. 

It is, however, possible, by inverting the line shape, to reconstruct the repulsive 
part of the dissociative potential in a coordinate range, over which the vibrational 
wavefunctions of the bound state are non negligible [36| . 

As only few vibrational levels of the bound state are thermally excited, the choice 
of the analytical form of the bound state potential is not too critical and we assume 
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that the bound state is correctly described by the theoretical literature potential 
(3)0+ 15]. We also assume that the repulsive part of the potential of the dissociative 
(1)0+ state is described by an inverse power law potential of the type 



V r = A+^ {x = R/R e ,) 



(7) 



in which A and Vq are adjustable parameters. For each pair (A, Vq) the scattering 
eigenfunctions are to be computed anew in order to calculate the Franck-Condon 
factors. The best agreement with the observed spectrum is obtained by setting 
A = (5315 ± 32) cm -1 and Vo — (760 ± 16) cm -1 . The spectrum computed by 
using these values of the parameters is shown as a thick solid line in Fig. [5] and 
is in nearly perfect agreement with the experimental one. The blue asymmetry of the 
experimental spectrum is also well reproduced and is due to the contributions of the 
vibrational states with v' > 0. 
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Figure 3. Comparison of the potential energy curves for the dissociative state 
O 4 ^ . Thin solid line: literature potential [15| . Thick solid line: best fit potential. 
The upper bound state potential (3)0j is reported along with two vibrational 
eigenfunctions for v' = and 8. 



In Fig. [3]we compare the theoretical and best fit potentials. We also plot the potential 
energy curve of the upper bound state along with the vibrational eigenfunctions for 
v' = and 8, respectively, in order to show over what coordinate range the comparison 
between the two determinations of the dissociative potential is meaningful. 

The difference between the theoretical and best fit potentials (approximately 
880 cm -1 for R = R e >) is fairly large. Probably, a theoretical determination of the 
electronic transition moments for these states or an improvement of the theoretical 
calculations of the potential energy curves might reduce the disagreement between the 
experimental and theoretical determination of the potential energy curves. 
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In any case, it has been demonstrated that the IR emission spectrum is due to a 
transition between two molecular states correlated with atomic limits in the 6p and 
6s manifold, respectively. 



3. High density results 

When the density N of the gas is increased, the increasingly stronger interaction of 
the Xe2 excimers with the environment affects the features of the emission spectrum. 
On one hand, as expected, the collision frequency increase reduces the lifetime of this 
species. Thus, the spectrum width broadens linearly with N. 




Figure 4. Density-dependent shift of the center of the excimer emission 
spectrum. Open dots: pure Xe. Closed dots: Ar(90%)-Xe(10%) mixture. Lines: 
prediction of the theoretical model. 



On the other hand, the center of the emission band shifts unexpectedly to lower wave 
numbers as N increases. A linear relationship of the type 

v m = 1>m,0 - AN (8) 

describes accurately the experimental data. 

This linear trend is observed both in pure Xe gas and in the mixture, although the 
slope has quite different values in the two different gases, and is obeyed for densities 
as high as TV 8 • 10 26 m~ 3 , i. e., for densities nearly 40 times larger than the density 
of the ideal gas at T = 273 K and P — 0.1 MPa. 

In Fig. [?] the central wave number of the spectrum, v m , is shown as a function 
of N. In the limit N — > 0, v m has the same value, the Xe2 excimer being produced 
cither in the pure gas or in the mixture. 

By contrast, when the Xe2 excimer is produced in the pure gas, the density 
dependence of v m is much stronger than in the case in which Xe2 is produced in the 
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mixture. If the data are fit to straight lines, the following values for the slope are 
obtained: A = (2.08 ± 0.04) • 10~ 22 to 2 in pure Xe and A = (0.63 ± 0.015) • 10~ 22 
m~ 2 in the mixture. We have to stress the fact, however, that the lines in Fig. 2] 
are not such linear fits but represent the predictions of a simple theoretical model 
that we have developed in order to understand the physics of the phenomenon we are 
investigating [291 ] . 

According to the suggestions of Mulliken [H] , an excited state of a diatomic rare 
gas molecule can be envisaged as an ionic core, e. g., Xe^~, plus a Rydberg electron in 
an extended orbital. We assume that the electron is so largely delocalized that many 
atoms of the host gas are encompassed within the large electron orbit. These atoms 
are polarized by the strong field of the ionic core. If the density is large enough, we 
can assume that the dielectric screening exerted by the polarized atoms can be treated 
by introducing the density-dependent dielectric constant of the gas, K(N). 

The interaction potentials are of electrostatic nature and are therefore reduced 
by a factor l/K 2 with respect to the dilute gas condition. At the densities of the 
experiment, the dielectric constant can be Taylor-expanded to first order as a function 
of N : K ks 1 + Na/eo, where a is the atomic polarizability of the atoms of the host 
gas and eo it the vacuum permittivity. 

The energy difference between the upper bound potential V + and the lower 
dissociative potential V Q + is reduced by the same factor. The average energy difference 
that approximately corresponds to the central wave number of the observed emission 
spectrum is thus given by 



hcK 2 (N) --'^y eo 

in which u m fi = (V + —V + ) / he and (...) is a suitable weighted average of the potentials 
over the coordinates of the vibrational states of the bound potentials. v m ,o is the 
zero-density limit of the central wave number of the emission spectrum. A linear 
dependence of v m is thus obtained. 

Unfortunately, if only this solvation effect were accounted for, the slope v m fi a l e ^ 
would amount to only half of the observed value. The missing contribution arises 
from quantum many-body effects in the interaction of the delocalized electron with 
the atoms of the surrounding gas. In fact, the wavelength of the loosely bound Rydberg 
electron in the excited dimer is very large and spans over many surroundings atoms. 
Such a state can exist provided that the electron is weakly scattered off the gas atoms, 
i. e., if the electron mean free path is much larger than the radius of its orbit, as it 
occurs in the present case fijl ]. 

In this situation, the Rydberg electron can be considered as nearly free and its 
simultaneous interaction with many atoms of the host gas leads to a density-dependent 
shift of its energy given by [43j 

AE(N) = ^-Na (10) 

TO 

where a is the electron-atom scattering length and to is the electron mass. 

This contribution has to be added to the energy difference between the initial 
and final state of the molecule, eventually leading to the following expression for the 
energy of the emitted IR photon in the center of the spectrum 

/ 2a \ h ( 2a ha\ , . 

v m = Vmfl 1 N H Na = u m - u m N (11) 

V e / toc V eo mcJ 
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The present model predicts a linear density-dependence of the wave number 
corresponding to the maximum of the spectrum, as actually observed. Moreover, 
the slope can be calculated if the gas parameters a and a are known. In the case of 
the heavier rare gases, the scattering length is negative and its contribution effectively 
adds to the solvation contribution. 

In the pure gas, the choice of the gas parameters is obvious. In the case of the 
mixture, the low concentration of Xe leads to the assumption that the excimer is 
actually surrounded by Ar atoms so that the parameters of Ar have to be chosen. 



In the case of the pure Xe gas, a = 4.45 • 10~ 45 F-m [4J] and a = -3.09 A |45|], 
yielding (2£ mi oa/eo — ha/mc) = 1.98T0 -22 m 2 , to be compared with the experimental 
value (2.08 ± 0.04) • 10~ 22 m 2 . In the cas of the mixture, the Ar parameters are 
a = 1.827 • 10~ 40 F-m [44] and a = -0.86 A yielding a slope value 0.66 • 10~ 22 
m 2 , to be compared with the experimental value (0.63 ± 0.015) • 10~ 22 m~ 2 . 

The theoretical prediction agrees very well with the experiment. In Fig fj] the 
straight lines actually represent the theoretical prediction and are not a mere linear 
fit to the data. 

It is interesting to note that the dielectric screening effect always produces a red- 
shift of the excimer emission spectrum because it reduces the strength of the Coulomb 
interaction between the Rydberg electron and the ion core. On the contrary, the 
multiple scattering effect gives a contribution to the energy of the emitted photon 
whose sign depends on the nature of the (nearly free) electron-atom interaction. If 
this is attractive, as in the case of the heavier rare gases, the electron-atom scattering 
length is negative, whereas it is positive in the case of Ne and He. Thus, roughly 
speaking, the density-dependent shift of the emission spectrum might be "fine tuned" 
from red- to (even) blue-shift by suitably choosing the buffer gas in which Xe excimers 
are produced. Further investigations in different mixtures should help confirming this 
model. 



4. Conclusions 

The study of the infrared emission spectrum of Xe excimers in pure Xe gas or in 
an Ar-Xe mixture reveals interesting features in both the low- and high-density 
limits. The analysis of the spectrum at low density has allowed the assignment of 
the molecular states involved in the transition. This was never done before because 
no IR molecular spectra of Xe in the present range were observed before. Moreover, it 
has been possible to give an experimental determination of the potential energy curve 
of a gerade molecular state related to the 6s manifold. 

On the other hand, it has been possible to see how cooperative- and many-body 
effects continuously develop by increasing the density up to quite high values. The 
change of the excimer emission spectrum with the gas density reveals that (at least) 
two kinds of mechanisms, one classical and one quantum, affect the interaction of the 
excimer with its environment. In our opinion, it is worth noting that the quantum 
contribution to the self-energy of the nearly free Rydberg electron links a typical issue 
of molecular physics such as the emission spectrum of a decaying diatomic molecule 
to the (apparently) uncorrelated realm of electron transport in dense gases. In fact, 
also in that case the density dependent shift of the electron energy due to multiple 



scattering deeply affects electron transport [46|; |47 
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